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Abstract 

The Richtmyer-Meshkov instability (RMI) refers to the baroclinic 
generation of vorticity at a perturbed density interface when impacted 
by a shock wave. It is often thought of as the impulsive limit of the 
Ray leigh- Taylor instability (RTI). The fluid dynamics video |"large-| 
eddy simulat ions (LES) of RMI in a converging geometry" | shows the 
mixing of materials resulting from the interaction of an imploding 
cylindrical shock wave with a concentric perturbed interface that sepa- 
rates outside light gas from heavy gas (initially 5 times denser) inside a 
wedge. At the initial impact, the incident shock Mach number is either 
1.3 or 2.0. The present canonical simulations support recent interests 
on compressible turbulent mixing in converging geometries relevant to 
both inertial confinement fusion and core-collapse supernovae dynam- 
ics. 

1 Flow configuration 

A converging cylindrical shock initially impacts a perturbed cylindrically- 
shaped density interface that separates (light) air from (heavy) SFq (density 
ratio of 5), both at rest at a radial position Rq = lm. Pressure and tem- 
perature are initially continuous across the interface, and set to 23kPa and 
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286K, similarly to the experiments of Vetter & Sturtevant in plane geome- 
try [I]. Two incident shock Mach numbers Mj = 1.3 and 2.0, defined just 
before shock-interface refraction, are tested. 

The flow is confined to a 90° wedge. We impose periodic boundary 
conditions in the z-direction of the cylinder axis (axial extent L z = lm). 
Zero-gradient boundary conditions are employed over an outer cylindrical 
boundary. It is assumed that shock wave-boundary layer interaction does 
not play a dominant role in the growth of the mixing layer and slip boundary 
conditions can be applied at the horizontal and vertical reflecting walls. An 
inner cylindrical wall is also used to regularize the apex. 

Initial irregularities in the density interface form the misalignment be- 
tween density and pressure gradients required to initiate RMI. A second 
RMI occurs after the initial shock has converged down the wedge, reflected 
off the axis and reshocks the distorted interface. Reshock interactions of de- 
creasing intensity follow successively. Due to the converging geometry, the 
accelerated or decelerated motion of the interface also generates RTI. As 
the fluids start to interpenetrate, secondary Kelvin-Helmholtz instabilities 
develop along the sides of the fingering structures. 

2 Numerical approach 

The reshock produces a large dynamical range of turbulent scales, requir- 
ing the techniques of LES. We employed the stretched- vortex subgrid-scale 
model of turbulent and scalar transport based on an explicit structural mod- 
eling of small-scale dynamics, as developed by Pullin and co-workers. The 
imploding nature of the flow is particularly suitable for the use of adaptive 
mesh refinement (AMR) provided by the parallel block-structured AMR 
framework AMROC. The Favre-filtered Navier-Stokes equations are solved 
on each Cartesian uniform subgrid of the mesh hierarchy. A weighted, essen- 
tially non-oscillatory scheme is used to capture discontinuities but reverts to 
a low- numerical dissipation, explicit, tuned center-difference stencil in the 
smooth or turbulent flow regions, optimal for the functioning of our explicit 
LES (See [2] for more details). 

The (x, y, z)-domain is discretized with 95 x 95 x 64 cubic cells on the 
base grid with three additional levels of refinement based on the local density 
gradient, reducing the computational expenses compared to the equivalent 
finest unigrid 760 x 760 x 512 problem. The initial perturbation wavelength 
is 100 times larger than the finest grid size, which is itself about 100 times 
larger than the estimated smallest physical scale in the flow - the Kolmogorov 
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scale - where the viscous dissipation of energy occurs. 



3 Video 

The video "LES of RMI in a converging geometry" shows four animations 



for each incident shock Mach number: 

1. magnitude of Schlieren density contours in perspective view; 

2. iso-surfaces of heavy- fluid mass fraction Ys~p 6 in perspective view (be- 
tween iso-surfaces = 0.1 and 0.9); 

3. levels of Schlieren density and heavy- fluid mass fraction in azimuthal 
and axial plane slices through the mixing region; 

4. AMR levels of refinement in perspective view. 

For the case Mj = 1.3 (resp. 2.0), each animation shows the flow evolution 
for about 25 xl0~ 3 s (resp. 15 xl0~ 3 s), and is played at a speed about 
800 (resp. 1,100) times slower than in experimental conditions. Data are 
sampled about every 0.5xl0~ 3 s. 

The first animation successively shows the location of the incident shock 
before impacting the initial perturbed interface; growing mushroom-like 
structures characteristic of baroclinic instabilities as the interface acceler- 
ates towards the axis of the cylinder, following the transmitted imploding 
shock; the reflection of the imploding wave off the axis and the first reshock 
interaction with the interface; the subsequent interpenetration of fluids and 
the development of a wide range of scales. The second movie portrays the 
mixing of light fluid (green) with heavy fluid (red) by displaying isosur- 
faces of heavy-fluid mass fraction. Both Schlieren density and heavy-fluid 
mass fraction representations are displayed together in the third sequence, 
where azimuthal and axial plane slices through the mixing region are set 
side by side, showing the contrast between azimuthal and axial modes. The 
fourth animation indicates the dynamics of each level of grid-refinement and 
suggests the computational savings achieved thanks to the AMR strategy. 
Strong incident shock strengths (e.g. Mj = 2.0) focus and compress the 
flow further down towards the axis, reshock interactions occur closer to the 
axis with greater intensity, leading to a stronger turbulent mixing. 
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